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INTRODUCTION 

The dynamic behavior of a helicopter rotor blade can 
be simulated using a P AC E® TR-48 Desk Top Analog 
Computer as the model. The computer with elec- 
tronic switching, solid-state components and high- 
speed repetitive operation solves the integrodiffer- 
ential and other equations describing the system. 
It then provides high-speed read-out of the results. 
These Notes provide a brief analysis of the system 
dynamics involved and a description of the com- 
puter model performing the simulation. 

SYSTEM ANALYSIS 

The primary result of the helicopter's main rotor 
action is a vertical lifting force. However, the ro- 
tor's motion also develops forces required for hori- 
zontal propulsion, aircraft stabilization and control. 
All of these forces result from the tilting of the rotor 
head or the entire aircraft, and from the 1 'flapping" 
of the rotor blades. The blades are attached to the 
rotor drive shaft by hinges that permit them to flap 
in the direction normal to the plane of rotation. 

When the blades are stationary, they droop down 
against stops. When in motion they turn upward, 
and are kept at a level close to the plane of rota- 
tion by a large centrifugal force. The forces acting 
on the rotor blades are indicated in figure 1. 
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Figure 1: Side View Of Helicopter Rotor Blade 



The analytic procedure used provides a set of 
equations that describe the dynamic behavior of a 
flapping rotor blade element. These equations de- 
fine parameters for the system force balance equa- 
tion. The latter, obtained by summing moments 
about the flapping axis, is the integrodifferential 
equation, K 

I (Sin 0Cos/3u +^f ) = 1/2 P C E |C L U rdr (1) 

o 

in which, 

C^ = lift coefficient, a non- analytic function of the 
angle of attack 

U =velocity of the blade element (in ft./sec.) 

Ce =mean cord in circle of rotation (in ft.) 

p =air density (in slugs/cu. ft.) 

I =blade moment of inertia (in slug-sq. ft.) 

u = rotor angular velocity (in radians/sec.) 

0 =blade flap angle (in radians) 

r =radius of blade element (in ft.) 

The angle of attack (a) is equal to the sum of the 
rotor pitch angle and the angle by which the hori- 
zontal (tangential) velocity is displaced from the 
true air speed. The rotor pitch angle is a function 
of the 4 * stick" position described by: 

A^g =the basic rotor pitch controlled by "con- 
trol stick" 

A^g=the pitch generated primarily by lateral 
movement of the "cyclic stick" 

B LS =the pitch generated primarily by longi- 
tudinal movement of the "cyclic stick". 

COMPUTER SIMULATION 

The details of blade motion are important in analyz- 
ing hinge-transmitted, shear stresses and other 
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causes of aircraft vibration. In addition, a knowledge 
of the blade dynamics is required in investigating 
stall as a function of ' 1 stick' ' position and angle of 
attack. The computer simulation of blade dynamics 
produces contour plots showing the angle of attack 
for various i ' stick" positions. 

A two- speed computation technique is used to solve 
the system equations. Continuous real-time compu- 
tation provides solutions to equations that are 
dependent upon azimuth angle (cot). These solutions 
in turn provide the inputs for computer circuits 
that operate with rotor blade length (X) as the in- 
dependent variable. High-speed repetitive opera- 
tion of these circuits provides the integral input to 
the force balance equation. 
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The TR-48 Computer model uses 46 amplifiers in 
addition to high-speed electronic- switches, com- 
parators and multipliers. A variable diode function 
generator provides the lift coefficient vs. angle 
of attack function. 

RESULTS 

The simulation produced a set of contour plots 
shown in figures 2, 3 and 4. They show angle of 
attack (a) for the various ' ' stick' ' settings de- 
scribed by A A LS and B LS . The plots indicate 
that the helicopter would stall for control stick 
settings (A cs ) of 10° and 12.5°. For these settings, 
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Figure 3: Stall "Stick" Settings; A cs = 10°, A LS = 1.25°, 
B LS = 6.25° 



the desired initial angle of attack (a Q ) exceeds 12°, 
a critical value representing the start of the stall 
region in the C vs. a function. 

This application demonstrates a practical technique 
for solving integrodifferential equations with the 
desk- top analog computer and for effectively read- 
ing-out the high-speed results obtained. 
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Figure 4: Stall "Stick" Settings; A cs = 12.5°, A L$ 
B LS = 6.25° 
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Figure 2: Non-Stall "Stick" Settings; A c $ = 5°, A|_$ = 125°, 
B L S=6.25° 



Contour Plots Showing Angle of Attack (a) vs. Azi- 
muth Angle (^) For Various "Stick" Settings. 
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